Landau-Lifshitz theory of the spin Seebeck effect 
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The spin Seebeck effect, i.e., thermal spin angular momentum transfer between a paramagnetic 
metal and ferromagnetic insulator, is studied theoretically based on the stochastic Landau-Lifshitz- 
Gilbert (LLG) phenomenology. Magnons in the ferromagnet establish a nonequilibrium steady state 
by equilibrating with phonons via bulk Gilbert damping and electrons in the paramagnet via spin 
pumping, according to the fluctuation-dissipation theorem. Subthermal magnons and the associated 
spin currents are treated classically, while the appropriate quantum crossover is imposed on high- 
frequency magnetic fluctuations. We identify several lengthscales in the ferromagnet, which govern 
qualitative changes in the dependence of the spin Seebeck coefficient on the magnetic film thickness. 

PACS numbers: 85.75.-d,72.25.Mk,73.50.Lw,71.36.-(-c 



I. INTRODUCTION 

Over the past three decades, spintronics has evolved 
from a focus on equilibrium phenomena in magnetic het- 
erostructures, such as giant magnetoresistance^ and in- 
terlayer exchange interactions ,'2' to dynamic processes, 
such as spin-transfer torque^ and spin pumpingjI^Sland, 
more recently, nonequilibrium thermodynamics, heralded 
by the spin Seebeck effecpE! and thermally-induced mo- 
tion of domain walls.'^'^ From a practical standpoint, 
magnetic nanostructures are useful for field sensing and 
nonvolatile information storagejl^ where magnetoresis- 
tance is paramount for the readout, while current- 
induced spin torques are useful for fast and scalable bit 
switching.E3 

One rapidly-developing avenue of research concerns 
out-of-equilibrium spin phenomena in insulating systems, 
where spin is carried by collective excitations, such as 
spin waves (magnons), rather than electronic quasipar- 
ticles. To this end, spin waves in the ferrimagnetic in- 
sulator yttrium iron garnet (YIG) appear particularly 
promising as they suffer from a remarkably low Gilbert 
damping (at microwave frequencies), a ~ 10""*, and 
the host material has Curie temperature of ~ 500 K, 
thus remaining magnetic at room temperature.121 Spin 
waves in YIG have recently been shown to undergo room- 
temperature Bose-Einstein condensation under nonlinear 
microwave pumping,i^ exhibit large spin pumping into 
adjacent conductors j'iS'i^ manifest the longitudinal spin 
Seebeck effect ,02 and efficiently move domain walls un- 
der small thermal gradients.^*' These phenomena hold 
promise for integrated circuits based on nonvolatile mag- 
netic elements^^ with essentially no Ohmic losses and 
thus very low dissipation. 

Furthermore, thermal control of magnetic dynamics 
and spin currents^ provides an attractive alternative to 
voltage control, especially since magnons, which are neu- 
tral objects, can respond more directly to temperature 
gradients. The spin Seebeck effect, i.e., the generation 
of thermal spin current between magnetic insulators and 
normal metals, is the basic phenomenon of central in- 
terest in this context. The purpose of this paper is to 



develop a rigorous phenomenological approach to this 
problem, based on the Landau-Lifshitz-Gilbert (LLG) 
phenomenology of ferromagnetic dynamics. ^^ Conceptu- 
ally, we are departing from the spin-pumping perspec- 
tive on the interaction between electrons and magnons at 
ferromagnetic-insulator I normal- metal interfaces put for- 
ward in Ref. |8j In this paper, however, we develop a sys- 
tematic Landau-Lifshitz treatment for thermal magnons, 
rather than monodomain dynamics central to Ref. 8., In 



Sec. VIII we comment on how the theory could be ex- 
panded to account for magnon and phonon kinetics when 
the standard LLG phenomenology fails. 



II. FERROMAGNETIC BULK DYNAMICS 

In the ferromagnetic bulk, the equation of motion is 
given by the stochastic LLG equatiorPil 

dtm = jm X (Hoff + h/) + am x dtui , (1) 

where m = M./Ms is the unit-vector magnetization di- 
rection {Ms = |M| being the saturated magnetization 
magnitude), 7 (minus) the gyromagnetic ratio (7 > 
for free electrons), a dimensionless Gilbert damping con- 
stant. 



H 



off' 



-^M^ = HaZ + A^W^m + Hj. 



(2) 



the effective field (consisting of applied field Ha in the 
z direction, exchange field ex A^, and relativistic cor- 
rections Hr that include dipolar interactions and crys- 
talline anisotropics), and h; random Langevin field with 
correlatop^S 



{hi4r,t)hi,j{r',t')) 



2a 



kBT{r)S,,S{r~r')S{t-t'), 

(3) 

in accordance with the fluctuation-dissipation theorem. 
We are interested at intermediate temperatures: much 
lower than the Curie temperature, such that the Landau- 
Lifshitz phenomenology based on the directional magne- 
tization dynamics [SO (3) nonlinear a model] is appro- 
priate, while not too low such that the classical theory 
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FIG. 1. Schematic of an Ni|F|N2 sandwich structure studied 
in this paper. The normal-metal layer Ni is treated as a poor 
spin sink, which blocks spin current, jsi ~ 0. The normal- 
metal layer N2, on the other hand, is a perfect spin sink, thus 
establishing a thermal contact between its itinerant electrons 
and magnons in the ferromagnetic insulator (F), which results 
in spin current J32 = js- We assume the phonons in the F layer 
follow a linear temperature profile from Ti at x — (Ni|F 
interface) to T2 at x = d (F|N2 interface), corresponding to 
electron temperatures in Ni and N2, respectively. 



can be used as a starting point. We will, furthermore, 
neglect H,. in Eq. (pi), for simplicity, which would be jus- 
tified when ksT ^ h^Ms- The Langevin correlator (pi) 
is white at frequencies w ^ ksT/h, corresponding to 
classical behavior. In Sec. |VI[ we will adapt our theory 
to account for quantum fluctuations at oj > fc^T/fi, by 
matching with the fully quantum treatment of Ref. 1231 
The temperature T(r) in Eq. (pi) is taken to correspond 
to the x-dependent phonon temperature. 



T{x) = T, + -{T2-T,), 



(4) 



assuming Gilbert damping stems from the local magnon- 
phonon scattering. (We will return to discuss these as- 
sumptions in Sec. VIII ) A schematic of our layered sys- 
tem is shown in Fig. Ill 

In order to simplify discussion of the spin transfer at 
the F|N interfaces, let us switch from the magnetization 
to the spin density: 



1. 



-M 



AI, 



sn . 
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where s = Mg/j is the saturated spin density and n 
— m its direction. The LLG equation then becomes 



(5) 



s(l -I- anx)dtn + n x {Hz + h) + i9Js 



0, 



where 



(6) 



(7) 



is,i = -An X diU 

is identified as the magnetic spin current and 

{h,{r,t)hj{r\t')) = 2askBTir)S,jS{r~r')5{t--t'), (8) 

where i and j stand for the Cartesian coordinates. Here, 
H = MgHa and A = M^A^. In equilibrium, n == — z, 
assuming the applied field Ha > 0. 



III. BOUNDARY CONDITIONS 

The boundary conditions for a ferromagnet sand- 
wiched between two normal metals need to be simi- 
larly constructed to account both for deterministicP and 
stochasticPS spin-transfer torques. We will start with 
the former and then include the latter according to the 
fluctuation-dissipation theorem. 

We assume the spin current is blocked by the Ni layer 
at a; = 0, due to its weak spin-relaxation rate: 



= (a; = 0) . 



(9) 



In other words, the spin current pumped across the F|Ni 
interface is balanced by an equal backflow.'^ For our pur- 
poses, Ni can thus be replaced by an insulator, as long 
as it makes a good thermal contact with phonons in the 
ferromagnet. A net spin current across the F|N2 inter- 
face, on the other hand, is allowed, if we treat N2 as a 
perfect spin sink.l^ 



Hq^^ dn 

^— n X -- 

47r dt 



{x — d) . 



(10) 



where g"^^ is the real part of the dimensionless interfacial 
spin- mixing conductance (per unit area);^- We neglect 
the imaginary part of g^"^, which is expected to have 
a negligible effect on our flndings in virtually all cases 
of interest. A specific realization for such an Ni|F|N2 
trilayer could be provided by the Cu|YIG|Pt combina- 
tion, where Cu (Pt) is a light (heavy) element with weak 
(strong) spin-orbit interaction. (We will generalize our 
flndings to arbitrary Ni|F|N2 trilayers, such as symmetric 
Pt|YIG|Pt type structures or general asymmetric struc- 
tures, in Sec. VII ) 



For the two Ni|F|N2 interfaces, we correspondingly 
have the following (deterministic) boundary conditions 
for magnetic dynamics (as T ^^ 0): 



d^n = , 

AdT-n H dtn 

47r 



x = 
, X — d 



(11) 



reflecting continuity of spin current, which is given by 
Eq. (pi) inside the ferromagnet and Eqs. (pi) and ( 10 ) in 
Ni and N2, respectively, across the corresponding inter- 
faces. Since spin pumping (pO|_affects magnetic dynamics 
similarly to Gilbert damping^ it is accompanied with a 
similar stochastic term.'23] The latter can be accounted 
for by modifying the boundary condition at x = d: 



Ad^n + ^9fn + h' = , 

47r 



(12) 



where 



{hi{p,t)hr{p\t')) = ^kBT,d,,S{p~p'nt-t') (13) 

Zn 

and p = (2/, z) is the two-dimensional position along the 
interface at x = d. The Langevin correlator strength is 



proportional to the electron temperature T2 at the F|N2 
interface, since the noise originates in the thermal fluc- 
tuations of electronic spin currents in N2 . 

The spin Seebeck effect is embodied in the thermal- 
averaged spin current flowing through the F|N2 
interfaced 

j,,^ = -An X a^n ^ n X ( -^dtn + h' j . (14) 

Utilizing axial symmetry of our system with respect to 
the z axis and linearizing small-angle dynamics, it is con- 
venient to switch to complex notation: n = n^ — iuy. 
Thermal spin-current density, (j^.x) = jsZ, can thus be 
written as 



js = Aim {n* dxn)\^^^ . 



(15) 



Exploiting, furthermore, translational invariance in the 
yz plane, we find in the steady state: 



is ^ Alva 
where 



(27r)3 (27r)3,5(q-q')5(w-w') 



n(q,w)= f d'^pdte''''^'-'if'>n{p,d,t) 



(16) 



(17) 



is the Fourier transform over p and time t. Similarly 
transforming Langevin correlators, Eqs. ([8]) and (13 1, we 
have: 

{h{x, q, uj)*h{x', q', w')) =4(27r)3asfcsr(x) 

x6ix-x')S{(i-q')S{Lu-uj'), (18) 

for the bulk and 

{h\q,uj)*h\q',Lu')) :^ 4{2Trf a'skBT25{q- q')S{LO - Lo') 

for the F|N2 interface, defining 



(19) 



Ans 



(20) 



which has dimensions of length, a'/d is the enhanced 
Gilbert damping for a monodomain precession of the fer- 
romagnetic film.— 



IV. SPIN SEEBECK COEFFICIENT 

We now have all the necessary ingredients in order 
to evaluate the (longitudinal) spin Seebeck coefficient 
(which has units of inverse length squared)p51 



S: 



Js 



kB{Ti-T2) 



(21) 



of the Ni|F|N2 structure shown in Fig. [T| To simplify 
our subsequent analysis, let us optimize the notation, as 



follows. The stochastic LLG equation ^ in the film bulk 
is written as 



A{d^ — n'^)n{x, q, uj) — h{x, q, w) , 



(22) 



after linearizing transverse dynamics and Fourier trans- 
forming it in the yz plane and time. Here, 



H — {1 + ia)suj 
A 



(23) 



The stochastic boundary condition ed, x — d, Eq. (12), in 
this notation is 

A{dx~ K')n{x,q,Lu) ^ -h'{q,uj) {x ^ d) , (24) 

while dxTi = at X = 0. Here, 



ex SijJ 



(25) 



Eqs. (22)-(25l now form a closed system of inhomo- 



geneous linear differential equations, with source terms 
given by stochastic fields h and h' . These are straight- 
forward to solve for n using Green's functions. Substitut- 
ing the solution for n into Eq. (16 1, we find, after some 



algebra, the spin Seebeck coefficient (21) 



S 



aa' s^ 


1 d\ 

— 00 


2tt^ A^d j 




.d 
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/ dxx 




Jo 



dLOUJ 



cosh[K(a; — d)] 



K sinh(K(i) — k' cosh{Kd) 



(26) 



Integrating over the longitudinal coordinate x, this finally 
becomes 



S = 



Sn^A^d 



,., , diocu 

d q / : 2 

-00 \Ksinh{Kd) — K'cosh(Kd)| 



sin (Kid) sinh (K,.(i) 



(27) 



where k^ (k^) is the real (imaginary) part of k. Eq. (27) 



is our central results and the main departure point for 
the subsequent analysis. 

Let us, for convenience, define the following length- 
scales: 



is the magnetic exchange length, 

a' 



I' 



a 



(28) 



(29) 



is the spin-pumping length (i.e., the F thickness at which 
the monodomain Gilbert damping enhancement due to 
spin pumping^ equals the intrinsic damping). 



A 



hA 
sksT 



(30) 



is the thermal de Broglie wavelength (in the absence of 
applied field), where T is the ambient temperature, and 



1 = 



A 



a 



(31) 



is the decay length for thermal magnons in the bulk. In 
this notation, 



2 2 1 1 + ia hio 



(32) 



1.5 



i^ 



A_ 



1.0 



1.2 



V. QUASIPARTICLE APPROXIMATION 

In the following, we are primarily interested in the 
thickness, d, dependence of the spin Seebeck coefficient, 
S, assuming the the length-scale hierarchy A ^ /' ^ Z. 
In YIG, for example, taking 47rM, ~ 2 kG, A - 
1/2 X 10~^ erg/cm, and a ~ 10~^, we find the follow- 
ing lengths: (1) A < 1 nm, at room temperature, (2) 
I' ~ 100 nm, taking g'^^ ^ 10^^ cm^^ from Ref. [TBI and 
proportionately larger /' with g^^ ^ 5 x lO^'* cm~^ from 
Ref. [T7] {g^^ is very sensitive to the preparation and qual- 
ity of the YIG|metal interfaces), (3) I > 10 /im, and (4) 
^ ^ 10 nm at 1 kG (corresponding to typical magneto- 
static fields). 

We start by performing integration over frequency uj 
in Eq. (27) in the limit of low damping (both intrinsic 
and spin pumping). In this case, the integrand is peaked 
at sinh(K(i) w 0, corresponding to 



1.4 1.6 

oj/oJo{q) 



2.0 



FIG. 2. Schematic plot of the integrand of the spin Seebeck 
coefficient, Eq. (27 1, illustrating the first four quasiparticle 
resonances, n = 0,1,2,3, according to Eq. (33 1, for a fixed 
q. We set a = 10"'', a' = 10"^ nm, g^ -f f^ = 1 nm~^, 
and d = 10 nm in the main plot. The inset shows the essen- 



tially continuum spectrum when d 
parameters unmodified. 



100 ^m, keeping other 



which can be readily integrated over frequency: 



So 



An^d J_^ {a + oi' Id) [1 + (a + a'/rf)2] 
cJ Id p d^q 



1 + VId }_^ (27r) 



(35) 



where we have assumed small damping, a -\- cJ Id <^ 1. 
Similarly, we find for the n > subbandsP^ 



oj^{g) = -\g +^^ + ^ 



(33) 



Or, 



la' Id 



d^q 



1 + 2Vld J_^ (2^)2 • 



(36) 



The total Seebeck coefhcient in the quasiparticle approx- 
imation is thus 



5* — 5*0 + 2_^ ^n 



(37) 



r!>0 



These resonances are well separated when their width is 
much smaller than their spacing, allowing for a quasipar- 
ticle treatment of the energy integral. For the bulk damp- 
ing, this condition is a ^ l/d-^/l/A^ — 1/^^^ for thermal 
magnons. Additionally, the occupation of magnons is ex- 
ponentially suppressed when the temperature is smaller 

than the gap, i.e., UbTIU < wq = His. Therefore, we are The wave vector q in integrals ([35|) and (|36| is bounded 
interested in the opposite regime, when the thermal de 
Broglie wavelength is smaller than the magnetic exchange 
length, A < ^. Thus in the regime d -^ I, these reso- 
nances are well-defined quasiparticle peaks corresponding 
to monodomain precession {n = 0) and standing waves 
(n > 0) along the longitudinal direction. See Fig. [2] This 
allows us to evaluate the spin Seebeck coefficient by sum- 
ming the contributions from individual magnon modes, 
when d <S^ I. 

Expanding around wq, the contribution from the lowest 
energy resonance is 



by requiring that huJn{q) < ksT, at which point our clas- 
sical treatment breaks down, as discussed in the next sec- 
tion. The spurious ultraviolet divergence is eliminated 
by cutting off at a total three-dimensional wave num- 
ber, Qc = •\/l/A2 — 1/^2^ corresponding to energy ksT. 
While a more careful quantum treatment for large wave 
numbers is constructed in the next section, we expect a 
crude cutoff to adequately capture the behavior of the 
Seebeck coefficient, up to an overall factor of order one. 
Allow us to momentarily focus on the regime when d ^ 
A, where the quasiparticle peaks are dense (recovering 
three-dimensional behavior), so that 



aa 



d^q 



dujuj 



S = 

Air^d J_^ " ^ i_oc {uj - ujo)'^ + {a + a'Idf a;2 ' 



(34) 



5" 



n>0 



bn 



2a' 



dh 



1 + 21' Id J (27r)3 l + 21'ldi'iT'^ 



(38) 



and, therefore, 




d<r 



d>r 



(39) 



That is, when the thickness of the ferromagnet is much 
smaller or larger than spin-pumping length, the spin cur- 
rent scales linearly with d or is d independent, respec- 
tively. 

If the quasiparticle peaks are not dense, d ~ A, finite- 
size effects are important, reflecting individual magnon 
subbands. In the extreme low-temperature case when 
d <^ X, only monodomain precession along the longi- 
tudinal direction contributes to the Seebeck coefficient: 
S Ki Sq. If the transverse dimensions are also much 
smaller than A, the full volume of the (nano)magnct un- 
dergoes stochastic monodomain precession, and 



S 



a 



1 



Vl + l'/d 



(40) 



where we have retained only one mode associated with 
the transverse momentum in Eq. ( 35 1 . V here is the vol- 



ume of the F layer. This coincides with the spin Seebeck 
coefficient for a monodomain obtained in Ref. Bl [defining 
(Ti + T2)/2 -^ Tp, Ta -^ Tn, and a'/d -^ a', to match 
their notation]. 

Finally, for largest thicknesses d^ I, the quasiparticles 
are no longer well-defined (see inset of Fig. [2]) and the 
above analysis cannot be applied. Because the thickness 
is beyond the magnon propagation length, only magnons 
within a distance I from the FjNa interface contribute to 
the spin current, which should, therefore, be independent 
of thickness, d, for a fixed thermal gradient, (Ti — T2)/d. 
Since, in this regime, the magnon propagation length is 
the largest length scale in the problem, we can send d — >■ 
cxD in Eq. ( 27 ) , which gives 



S^ 



ilT^A^d 



d^q 



dujuj 



oo ^^ |/^ ^1 



(41) 



The integrand, which can be evaluated numerically, is 
independent of thickness, and, therefore, S oc 1/d, as 
expected. 



To this end, we use a quantum-mechanical resullP^ for 
the thermal spin current, which is exact for a tunneling 
spin-exchange Hamiltonian at an F|N interface!^ 



js = -AkeSTa' / deeD{e)P 



,5riBE(/3e) 
9/3 



(42) 



where -D(e) is the magnon density of states, nBE{x) = 
(e^ — 1)^^ is the Bose-Einstein distribution, /? = l/ksT, 
and 5T is the temperature drop across the F|N inter- 
face (assuming magnons are equilibrated to a uniform 
temperature T, such that A ^ d). This limit can be di- 
rectly compared to our Eqs. ( plj ) and (27 1, by first send- 
ing a' -^ in the integrand (thus reproducing the weak 
F|N contact and allowing for the magnons in F to equili- 
brate with phonons) and then sending a — > [such that 
the magnon spectral properties are unaffected by Gilbert 
damping, as assumed in the derivation of Eq. (42)]. The 
magnons are correspondingly equilibrated to the average 
phonon temperature (Ti -I- T'2)/2, such that we identify 
ST = (Ti — T2)/2, in the present notation. According to 
Eq. ( [38| , our semiclassical spin current becomes in this 
limit 



js-^2kB{Ti~T2)a' 



d\ 

(27r)3 



(43) 



This agrees with Eq. (H2| in the limit e <^ ksT , where 



'^«Di.)fi'.^^^iy.M^,.l^ 



dp 



(44) 

We conclude that the classical-to-quantum crossover 
could be accounted for by inserting the factor 



e/3 



2 9nBE(/3e) 



dp 



(5huj/2 



smh{l3hijj/2) 



= FiPhuj) (45) 



in the energy integrand of Eq. ( |27[ ) , which effectively cuts 
off the contribution from magnons with energy e = hw ^ 
keT. 



VII. RESULTS 

We summarize the spin Seebeck coefficient dependence 
on the ferromagnetic layer thickness, when X <^ I' -^ I: 



VI. QUANTUM CROSSOVER 

Our classical Langevin theory needs to be appropri- 
ately modified when approaching magnon frequencies of 
hoj ~ ksT. On the one hand, this is an important limit, 
as the spin transport is dominated by thermal magnons 
in our model. On the other hand, the classical theory 
is inadequate for the treatment of quantum fluctuations 
that dominate at high (on the scale of the ambient tem- 
perature) frequencies. 



S{d) 





r a, 


d< A 


1 


d/l, 


A < d < r 


27rA2 ' 


I' /I, 


V ^d^l 




[I'/d, 


I < d 



(46) 



assuming A <C ^ (or else the magnon transport is ex- 
ponentially frozen out). To illustrate this dependence, 
we plot the spin Seebeck coefficient as a function of d 
in Fig. [3) using lengths characteristic of YIG, which are 
consistent with the above length-scale hierarchy. Notice 
that even though ^ determines the magnon energy gap 




FIG. 3. Plot of the spin Seebeck coefficient, Eq. (271, as a 



function of ferromagnet thickness d for A = 1 nm, ^ = 10 nm, 
I' = 100 nm, and I — 10 /jm. We use Eq. (37 1 (solid curve 
on the left) and Eq. (41 1 (solid line on the right), which are 
valid when d < I and d > I, respectively, when Ni is a poor 
spin sink and N2 a perfect spin sink. To account for the 
classical-to-quantum crossover, we have inserted factor ( |45[ ) 
in the integrands of Eqs. (351 and (361 [with oj — >■ u)„,{q)] and 
Eq. (41 1. The dotted curve shows the enhanced spin Seebeck 



coefficient when both normal layers are perfect spin sinks. 



and the associated Ginzburg-Landau correlation length 
in the classical theory, it does not govern any prominent 
crossover in the function S{d). 

We conclude that S{d) has nonmonotonic thickness de- 
pendence, with the maximum value 



^(r 



27rA3 ' 



(47) 



attained Sit I' < d < I, i.e., below the magnon decay 
length but above I', such that magnons equilibrate fully 
to the average phonon temperature f = {Ti + T2)/2 
(d > I') but still remain coherent on the scale oid {d < I). 
Amusingly, 5''™^'^-' is proportional to the spin-mixing con- 
ductance [see Eq. (20)] but is independent of the bulk 
Gilbert damping. According to Eq. (21), 5(™a-'^) deter- 



mines the largest spin current emitted thermally by a 
film of magnetic insulator, for any d, when subjected 
to a certain temperature difference (for example, in a 
wedged magnetic insulator coated by metallic contacts). 
If, on the other hand, a well-defined temperature gradi- 
ent is supplied (corresponding, for example, to a certain 
phonon-dominated heat-flux density), while thickness d 
is varied, the spin-current density js oc Sd increases with 
d saturating at d > I (the magnon decay length): 



J 



(max) 



I' 



fixed a^T 27rA2 



ked^T, 



(48) 



.(max) 



which corresponds to the bulk regime, js vanishes 

when a' — >■ (no spin pumping) or a — >■ 00 (no magnetic 
dynamics). 



It is interesting to ask how the above results would 
modify if both Ni and N2 in our model (see Fig. [I]) were 
perfect spin sinks. For an inversion-symmetric struc- 
ture (e.g., Pt|YIG|Pt), spin currents at the two inter- 
faces must be equal, jsi = js2 = js- When d ^ I, we 



should recover the bulk limit (41 1, since the magnons de- 



cay before traversing the full width of the film (and thus 
the spin current at one interface should not be sensitive 
to the boundary condition at the other). When d <S^ I, 
however. So and S'„ entering Eq. (37 1 need to be modi- 
fied. To that end, we notice that the factor {l + l' /d)~^ in 
Eq. (35 ) reflects the difference between Tm,o, the effective 



temperature of the magnons, and T2, the temperature of 
the electrons in N2P 



TmQ — T2 — 



aT + {a'/d)T2 _ Ti - T2 1 



a'/d 



l + l'/d 



(49) 



Similarly for the n > subbands, the factor {1 + 21' /d) ^ 
in Eq. ( 36 1 stems from the effective magnon-electron tem- 



perature difference across the F|N2 interface of 
^ Ti-Ta 1 



2 1 + 21' /d 



(50) 



In the case when both Ni and N2 are perfect spin 
sinks, the effective magnon-electron temperature differ- 
ence driving spin current is given simply by T — T2 — 
(Ti — T2)/2 for all subbands. We account for this in- 
creased thermal gradient by dropping the factor (1 -|- 
l'/d)~^ on the right-hand side of Eq. (35 1 and likewise 
(1 + 2l'/d)^^ in Eq. (36). The corresponding enhance- 



ment of the spin Seebeck coefficient reverses the trend in 
(461) at d < r 



Eq. 



S{d) 



to give 



1 



27rA2 



a'/d, 


d< A 


I' /I, 


A< d< / 


I'/d, 


I < d 



(51) 



which is now monotonically decreasing with d, as plotted 
by the dotted line in Fig. [3J 

When the structure Ni|F|N2 is not mirror symmetric 
(either because the spin-mixing conductances or the spin- 
sink characteristics are different), which we characterize 
by different a'l and a'2 spin-pumping parameters at the 
two interfaces, we can repeat the above analysis for d <C ^, 
finding 



T 



r?j.,0 



To 



and 



T„ 



T2 



aT+ia'Jd)Ti + {a'Jd)T2 

a + a'l/d + a'2/d 
T1-T2 l + 2l'Jd 
2 1 + (/'i + l'^)/d 



T1-T2 1 + 4/Vd 
2 \ + 2{l'^+l'^)/d' 



(52) 
(53) 

(54) 



where we defined the spin-pumping lengths 1',^ = a'^^/a as- 
sociated with the left, i = 1, and right, i = 2, interfaces. 



We thus generalize the spin Seebeck contributions (at the 
F|N2 interface) from different magnon subbands to 



So ? 
and 

Sn>a 



l + 2l[/d 



d 1 + {l[+l'^)/d J _^{27ry 



F[l3hiooiq)] (55) 



l+4l[/d 



d 1 + 2{l[ + V^)/d 



d\ 



F[phuj^{q)] (56) 



in heu of Eqs. (35 1 and \im, respectively. The asymptotic 





r 2/{a\-^+a'^-^)d, 


d< A 
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2/{l[^'+l',-')l, 


A < d < ?i 


27rA2 


lyi, 


l[,l'2^d^l ' 




[ 1'2/d, 


?< d 

(57) 



and crossover trends are now given by (assuming A <^ 



Sid) 



which coincides with Eq. ( |5l[ ) when a'l = a'2 but has an 
additional shoulder- like feature at d ^ (/'^ +12)/ 2 . This 
feature makes S{d) nonmonotonic when a'^ < a'2- 



VIII. DISCUSSION 

In order to compare our Seebeck coefficient to Ref. [S] 
[beyond the monodomain limit, Eq. (40 1, where we have 
full agreement] 



in the more interesting case of macro- 
scopic films, we consider the limit I' <^ d <^ I, where our 
result [see Eq. (46)], S ~ a'/2TrX^ should be compared 
with their S ^ a' /Vc, where Vc is the "magnetic coher- 
ence volume" (defined as a magnetic monodomain that 
has the same magnetic fluctuations as local fluctuations 
in the infinite bulk). The divergence of Vc oc ujq in the 
absence of applied field highlights the incompatibility of 
their findings with ours for films with dimensions larger 
than the thermal de Broglie wavelength A. 

We note that our theory provides the most minimal- 
istic (internally consistent) application of the LLG phe- 
nomenology to the problem of the spin Seebeck effect, yet 
disregards magnon-magnon interactions. These become 
important at high temperatures, especially approach- 
ing the Curie temperature. Magnon-phonon interactions 
are included only insofar as a contribution to the total 
Gilbert damping. Elastic magnon scattering on impuri- 
ties, which would manifest as an inhomogeneous broad- 
ening of ferromagnetic-resonance linewidth, may be an 
important impediment to the thermally-induced spin cur- 
rents in disordered films. The bulk limit of spin current, 
Eq. (48), is reduced by disorder, as well as the magnon 



decay length describing the crossover to the bulk regime. 
When the magnon mean free path I* is shorter than our 
Gilbert damping decay length I, in particular, we expect 
for the effective decay length to be leg ^ VU* and jf^^^' 
in Eq. ( 48 ) reduce by a factor of I /leg ~ \/l/l* (assuming 
that I' <^ Ics, such that our length-scale hierarchy is un- 
changed). The Seebeck coefficient behavior ( [46| [as well 
as Eq. (47)], however, remain essentially intact up to the 
thickness d ~ ^cff- 

Finally, we want to comment on a possibility of nonlo- 
cal magnetic relaxation. In this paper, we have assumed 
that Gilbert damping is a local and isotropic tensor. 
The locality would be a reasonable approximation if the 
damping bottleneck was due to some local dynamic de- 
fects. In the case of YIG ,123 which is known for its highly 
coherent elastic properties, nonlocality of the bulk mag- 
netic relaxation could significantly modify our findings. 
First of all, this could introduce new phonon-dependent 
Icngthscales into the problem, which would show in the 
S{d) dependence. Standard long- wavelength ferromag- 
netic resonance on thick films would reveal damping a 
that could be very different from that of short-wavelength 
thermal magnons relevant here. An effective damping 
parameter a of thermal magnons (which may itself be 
thickness dependent) would result in the bulk crossover 
thickness of I ~ A/a ^ I. When d<l,we may still invoke 
the first three regimes of our findings, Eq. ( 46 ) (with a 
and I corresponding to the thermal-magnon inverse qual- 
ity factor and decay length, respectively, due to magnon- 
phonon scattering), which should, furthermore, be indif- 
ferent to the fact that the local temperature, Eq. (|4|, is 
not well defined for highly-coherent phonons. The rea- 
son for this is that, in these regimes, when d is below the 
magnon decay length I, only the average phonon tem- 
perature T is relevant for our theory. Our bulk regime, 
S{d) ex d~^, when d is larger than the magnon decay 
length, would, however, have to be considerably revised 
in the case of nonlocal magnetic relaxation, calling for a 
rigorous quantum-kinetic theory.'^S' 

In summary, we have developed a minimal Landau- 
Lifshitz theory of the longitudinal spin Seebeck effect, 
which calls for its systematic experimental study for 
the temperature and film-thickness dependence, which, 
in turn, may necessitate a more systematic microscopic 
quantum-kinetic theory. 
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